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Abstract 

We study the annihilation process pp — > 7r°e + e~ within a Regge framework, as a means to provide contraints on timelike 
nucleon form factors. We present results for the e + e~ angular distributions and the differential cross sections in kinematics 
which will be accessible by PANDA@FAIR. To check the consistency of the model, we first test the approach on the process 
of real photon production, pp — > 7r°7, where data in the energy range of 2.911 GeV< yfs < 3.686 GeV exist. We find that a 
Regge pole model is able to well reproduce the available data. The analysis is then extended to a timelike virtual photon in 
the final state. 
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I. INTRODUCTION 



The structure of the nucleon in terms of the electromagnetic form factors can be probed by the electromagnetic 
interaction. In the spacelike region of negative momentum-transfer q 2 the electric and magnetic form factors, Ge and 
Gm, can be investigated in elastic electron-proton scattering. The spacelike form factors, which provide information 
on the electromagnetic spatial distributions in the nucleon, have been studied extensively, using both unpolarized 
measurements as well as polarization experiments, for reviews see e.g. Refs. [H-Q. The corresponding crossed processes 
pp — > e + e~ or e + e~ — > pp allow to access the form factors in the timelike region (q 2 > 0), starting from the threshold 
q 2 = 4m^r, where m^r is the nucleon mass. A few measurements of the latter processes, mainly of the total cross 
section, exist. Only two experiments, performed at BABAR Q and LEAR ||, attempted an individual extraction 
of the timelike electric and magnetic form factors. New measurements of the timelike form factors are planned at 
BES-III and PANDA@FAIR Q and will explore the at present still largely uncharted timelike region in much greater 
detail. However, the timelike region below the threshold < q 2 < 4m^, which is denoted as the unphysical region, is 
not accessible by annihilation processes as pp — > e + e" or e + e~ — > pp. Nevertheless, to obtain a consistent picture of 
the nucleon, knowledge of the electromagnetic form factors over the full range of q 2 is needed. 

Therefore, in Refs. 0, H[ the annihilation process, where in addition a neutral pion is produced, pp — > n°e + e~ has 
been studied. Since the outgoing pion takes a part of the energy of the process, the production of a lepton pair with 
an invariant mass below the (p + p)-annihilation threshold is possible and thus this reaction can be used to study 
the electromagnetic form factors in the unphysical region. Moreover this reaction provides the possibility to access 
the relative phases of Ge and Gm, in contrast to cross section measurements of pp — > e + e~ , where only the modulus 
of the form factors can be extracted. However, a study of the timelike form factors from the pp —> ir°e + e~ process 
requires a theory in order to extrapolate to the nucleon pole. 

In Ref. 0, the process has been calculated using a Born diagram model, in which the reaction is described through 
the exchange of a single nucleon, the Feynman diagrams arc shown in Fig. [TJ A formula for the differential cross 
section has been presented, where the dependence on the lepton kinematics has been removed by integration over the 
full phase space of the leptonic subprocess. The calculation has been extended in Q for an exclusive experimental 
setup including the investigation of polarization obscrvablcs. For a high invariant mass lepton pair the process has 
been studied in [lol. fll| using the concept of transition distribution amplitudes. 

In this paper, we study the process pp — > n°e + e~ using Regge theory [l2l. [ill], and describe the process at high 
energy within a Regge pole model. This phenomenological approach has been sucessfully applied to electroproduction 
and photoproduction of pions and kaons, see e.g. (Til . Il5l |. In particular, it has been widely a ppli ed in order to 
extract tt + and K + electromagnetic form factors from the 7r + and K + electroproduction process [l6l [ItJ • To check 
the consistency of this model, we first test the approach on the process of real photon production, pp — > 7r°7, which 
has been measured at Fermilab [l8| , and compare the results of our calculation with the existing data. We find that 
the Regge pole model is able to describe the results of the measurements, whereas the Born diagram model cannot 
reproduce the available data. 

After specifying the Regge model for pp — > ir°j, we extend the approach to the annihilation process pp — > ir°e + e~ 
and investigate the cross section with regard to the determination of the form factors in the unphysical region within 
the kinematic range, which will be accessible by PANDA@FAIR. Furthermore a general expression of the cross section 
including the dependence on the full lepton pair kinematics is presented. 

The paper is organized as follows: in Section [Til wc analyse the process pp — > ir°-f. We introduce the Regge pole 
model and compare our results to the existing data. In Section lllll the approach is applied to pp — > 7r°e + e~. The 
cross section in terms of the angular distribution of the produced lepton pair is given and the angular dependence 
of the cross section for several kinematic setups is presented within the Regge model. Wc summarize our findings in 
Section ITVl 



II. THE pp ->■ 7T°7 PROCESS WITHIN A REGGE FRAMEWORK 

Before describing the theoretical model for pp — > 7r°e + e~ process, for which at present no data exist, we firstly 
consider the reaction 

P(Pl)+p(P2)^^°(q^+ 1 (q), (1) 

which has been measured at Fermilab [l8l ]. and will test our model predictions with the results of this experiment. 
Data of the angular dependence of the differential cross section da/dcosO^ is available in the center-of-mass (cm.) 
energy range of 2.911 GcV< y/s < 3.686 GcV. 



2 




FIG. 1. Born model for pp — > tt e + e described by a single nucleon exchange in the it-channel and t-channel diagrams. 



The annihilation process can be described by two kinematic invariants, e.g. two of the Mandclstam variables, 
defined as: 

s = (p 1 +p 2 ) 2 1 u = (pa - qir) 2 ', t = (p l -q 1T ) 2 . (2) 

We start with the model considered in Ref. where the process pp — > n°e + e has been described within a Born 
diagram model through the exchange of a single nucleon. The two corresponding Fcynman diagrams, which are shown 
in Fig. [TJ arc given by a it-channel and t-channel nucleon exchange. The amplitude of the diagrams can be written as 

T " = i £ m ' -Wa-, u =i £ m ' «(Pi) r 7WJv(«) ( 7 P !, 9 I,2 +mJV ) T^Jvfe) u(pa), 

y i \ " m Af / ^ 



T t = 72 £ M ' ^o 7 . t =72 £ M ' ^JVJvCfe) ——2 r 7WA?(9) ufo), 

y i \ L " t jv / 

where r 7 jviv (r^Ar) is the parametrization of the ^NN* (irNN*) vertices and describes the leptonic subprocess 
7* — > e + e _ . 

Within this approach off-shell effects of the exchanged nucleons were neglected, hence the 7iViV*-vertices are 
parametrized by the on-shell proton electromagnetic form factors, in terms of the Dirac and Pauli form factors F\ 
and F 2 given by 



2m n 



(4) 

To describe the ttNN* vertex we will study both cases of pseudoscalar as well as pseudovector irNN coupling: 

g„NN{ml) (A75 + (1 - A )%7T'ys) . (5) 

where A = 1 (A = 0) leads to a vertex with a purely pseudoscalar (pseudovector) coupling, with the pion-nucleon 
constant g^NNim 2 ;) = 13.02. 

In case of real photon production pp — > 7r°7, the amplitudes within the Born diagram model are obtained as 



Tu = s tl - M* u =e M (g, A 7 ) • v {Pi)T^ NN — s T„ NN u{p 2 ), 

\ u N J 

t * kav * I \ \ -I \r f 1 ■ (Qtt -Pi) +m N \ 

% = e li - M^o t =e^(q, A 7 ) • v( Pl ) T^nn g u(p 2 ), 

\ 1 m N J 



(6) 



where e(q, A 7 ) is the real photon polarization vector and the form factors Fi(q 2 = 0) = 1 and F 2 (q 2 = 0) = up = 1.79 
have been used in the description of Tjnn- 

The differential cross section da/dt for the process pp — > ir°-f can be expressed as: 



da 



-(-g^)K^ (7) 



dt 16-7rs(s — 4m^) 
where W l ^o 1 is the hadronic tensor, given by 
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FIG. 2. Born model for p + p — } 7T° + 7 described by the exchange of A(1232). 



Using the Born diagram model for pp — > 7r°7, we are not able to reproduce the results of the E760 experiment. The 
obtained cross section da /d cos 6^ is about 4 to 5 orders of magnitude larger then the data, depending on the value 
of the cm. -energy y/s. Simple fixes by introducing strong suppressions through off-shell form factors do not lead to 
a correct energy dependence of the cross sections. One can assume, that this model is not suitable to describe this 
process and thus the process pp — > 7r°e + e _ as well. 

Therefore we consider a Regge pole description, which is based on the exchange of dominant baryon Regge trajec- 
tories in the w-channel and i-channel. This approach is valid in the kinematic region at forward and backward angles, 
s ^> |f| and s \u\. In the kincmatical region s ~ —t ~ — u the reaction has been investigated within the framework 
of generalized distribution amplitudes fl9l ]. The dominant trajectories for the process pp — > it°j are the nucleon (N) 
trajectory and A-trajectory associated with the A(1232)-resonance. The amplitude for Regge trajectory exchange 
can be obtained from the Born diagram by replacing the usual Feynman propagator of the single exchanged particle 
by the so-called Regge propagator, while leaving the Feynman structure, giving rise to the residue of the Regge pole, 
unchanged. 

In case of an exchanged nucleon, the pole-like Feynman propagators of the it-channel and t-channcl, given by 
l/(u— rn 2 N ) and l/(t— mfj), are replaced in the following way 



m T D N° SS °( U > S ) - r 



a N {u)-± -ll(o»(u)+i) 



N T[a N (u) + ^] sin7r (aAr(w) 5; 

(9) 

s ajv(t)-i — iir(a N (t)+±) 



T[a N (t) + ±] N smn(a N (t) + §) 



where the nucleon trajectory aiy is of the form ccn{u) = \ + a' N (u — m 2 N ) with a' N = 0.97 GeV~ 2 . 

Besides including the nucleon Regge propagators, we consider the exchange of the A-trajectory. Starting from the 
Feynman diagrams in a Born model, illustrated in Fig. [2j the amplitudes of the exchange of the A(1232)-resonancc 
can be expressed by 



(10) 



xW-sW 3^ 3m|-r r ^(P 2 ) 



M t°y,t =(-9»v) ■ v(pi)T% NA j——2 (7 • Pt + m A ) 

/ 1 laPt.fi - IPPt.a 2p t a p tJj \ g 

x|. 9q ,--7^ ^ 3^1-/^^)' 

with p u = p2 — q-K, and pt = p\ — q-n- T-yNA and F^jva are describing the 7./VA-vertices and irN A- vertices, respectively, 
and can be parametrized through 

F a ■ /2 3e(m A + mjy) 2s aflp a 

Kna = - ^-r^Q^PA^Tl (12) 



4 



where pa is the 4- momentum of the intermediate A-state. The operator TjJ is the isospin 1/2 — » 3/2 transition 
operator, f„ denotes the pion decay constant and Ua — 2.85 the 7riVA coupling constant. In Eq. (fT2"]) the 7^A 
vertex is parametrized by the form factor gu (q 2 ), the magnetic dipole form factor representing the strength of the 
magnetic dipole N — > A transition. The electric quadrupole and Coulomb quadrupole terms have been neglected in 
our calculation since their contributions have been found to be of order of a few % [201 ■ We will use as 7-/V A coupling 
strength g M (0) = 3.02. 

The Fcynman propagators in Eq. (|10p are then replaced by the Rcggc propagators: 



1 



L» A CgS °(u,s) = — = pr TT^ 



'A r [qa(u) + |] sin7r (oa(m) - |) ' 

(13) 

-jRegge/. * / e 



A \ b '°J-^l /.x I 11 ""A 



r[a A (*) + 3] sin7r(a A (i)- 1)' 



where the A-Regge trajectory is of the form q;a(w) = § + a'&(u — m A ), with a' A = 0.9 GeV -2 . 

We reggeize the amplitude of the process in the following way, which ensures gauge invariance of the Rcggc model 
amplitudes: 

M% tU = D^ c (u, s)(u m%) [M u + M t ], 

(14) 

M%, t = D^(t, s)(t - m%) [M u + M t ], 

and analogous expressions for the A-exchange diagrams. In the limits u —> m 2 N or t — > m 2 N , the results of the Born 
diagram model are recovered, since one obtains D^ C6SC (u, s) ■ (u—m 2 N ) — > 1 and D^ S6C (t, s) ■ (t—m 2 N ) — > 1 in this limits, 
respectively. Note that the Regge amplitude implies s^>|f|,s3>|u|, so that both forward and backward regions are 
kinematically separated. Only in this limit there is no double counting in this procedure. In the intermediate angular 
region the results of the Regge framework should only be taken as indication that one is outside the range of the 
validity of such an approach. 

The results of the differential cross section da /d cos 6„ including iV-trajectory exchange are shown in Fig. [3] for 
several cm. energies ^fs using two different pion-nucleon couplings, pseudoscalar and pseudovector, corresponding to 
A = 1 and A = in Eq. © , respectively. One notices that the angular dependence of the data is well reproduced and 
the results in the forward and backward regions are slightly smaller than the data. They differ by a factor of ~ 3. 
Using pseudovector couplings we obtain a cross section, which is smaller in magnitude and has a different angular 
behavior for larger values of | cos n | , compared to the case of pseudoscalar coupling. The deviation is decreasing for 
increasing cm. energies. Since the model with pseudoscalar coupling gives a better fit to the data, we will use this 
coupling in the following calculations. 

We next discuss the inclusion of the A-exchange Regge trajectories. As for the A we are further away from the pole 
position than in the nucleon case, the description of the residues of the Regge poles through their on-shell couplings 
can be expected to be modified. We allow for such a reduction of the coupling strengths of the A Regge pole residue, 
leading to the amplitude 

U^ = \Y,\MX+F-M^\\ (15) 

spins 



where the parameter J- will be obtained by a fit to the data. 

In Fig. |H results for do jd cos 6^ including A^-trajectory exchange and (N + A)-trajectories exchange are presented 
as well as results using the approach of Eq. (fT5"j) . Fitting the available data leads to T ps 0.5. We find that the Regge 
model including A-trajectory exchange describes the available data very well. The results give a better description 
compared to the case when only ./V-trajectories are taken into account. In particular for larger values of s, the results 
of the iV-trajectory lie below the data. The cross section including the reduction factor T of the A-polc residue is in 
very good agreement with the experiment, especially in the region s 3> \t\, \u\. 

In order to test the applicability of the model, we present in Fig. \E\ the cross section da/dt, divided by the expected 
s-dependence of the leading Regge trajectory, s 2q a(*)-2^ ag a f unc tion of — t for two different cm. energies. One 
notices, that for — t — > (da /dt) / s 2aA ^~ 2 shows a behavior, which is approximately independent of s, as expected 
from Regge theory, and approaches a constant value. The existing cross section data, as indicated by the data taken 
at t/s = 3.686 GeV in Fig. [5j have not yet reached the region of such small values of —t, where an extrapolation of 
(da I dt) I s 2aA ^" 2 could be performed by a constant. 
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FIG. 3. Differential cross section da jd cos 6-n of pp — > tt 7 for different cm. energies y^s including nucleon Regge trajectory 
exchange in the f-channel and it-channel: the blue (solid) curve corresponds to a pseudoscalar coupling of the ttNN vertex 
(A = 1 in Eq. (JSJ), the green (dashed curve) describes a pseudovector coupling (A = 0). The data are from Ref. [l8| . 



III. THE pp -> n°e+e- PROCESS WITHIN A REGGE FRAMEWORK 

After specifying the Regge pole model, we study the process 

P(Pi)+p(P2) -> 7T°(g OT ) + 7 *(g) -> tt ^) + e-(fti) + e+(fc 2 ), 



(16) 



and will start with a model independent analysis of the annihilation cross section. The cross section for the process 
is defined as 



da 



1 



d 3 q n 



d 3 h 



d 3 k 2 



~^(pi-P2) 2 -m% W^Ej \{2nf2klJ \{2n) 3 2k° 2 ) 
x (2n) 4 6( Pl +P2-q*-k 1 - k 2 )\T\ 2 , 
where the amplitude T can be decomposed into a hadronic and a leptonic contribution: 

E lE I ( M %* ■ e *nfa> A ^)) 72 ^Mki) v{k 2 ) 



\T\< 



X-y— 0,±1 spins 



(17) 



(18) 



A^o 7 » is the amplitude of the hadronic process pp — > 7r°7*. Both contributions of Eq. (JTSJ) 

1 2 1 2 

|(A4^ . ■ e*(q, A 7 )| and \e v {q, A 7 )?2(fci) e7 l/ v(fc 2 )| , 



(19) 



are invariant, thus we can choose any reference frame for our calculation. The advantage of such a separation is that 
one can easily calculate the hadronic and leptonic processes in two different reference frames. 
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FIG. 4. Differential cross section da jd cos 9^ of pp — > 7r°7 for different cm. energies yfs\ blue (dotted) curve: iV-trajectory 
contribution; purple (solid) curve: cross section including (N + A) Regge trajectory exchange; green (dashed) curve: (N + A) 
contribution including a reduction of the A pole residue (J- ~ 0.5) according to Eq. (|15[1 . The data are from Ref. [3]. 



We consider the hadronic process in the cm. -frame of the nucleon pair, where the momenta of the incoming proton 
and antiproton have opposite directions, illustrated in Fig. El The lcptonic subprocess 7* — > e + e~ is computed in the 
7*-rest frame, with the 4-momentum of the virtual photon given by 



q = (vV, 0,0,0) (20) 

and where the lepton momenta can be written as 

h =— — (1, sm0 e+e - cos$ e+e ~, sin6» e+e - sin $ e+e - , cos 6» e+e - ) , 

h (21) 

kt2 = — — (l, — sin# e + e - cos<f> e + e - , — sin^g+g- sin$ e + e - , — cos^g+g-). 
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0.1 1 

-t [GeV 2 ] 

FIG. 5. Cross section da/dt divided by s 2a& ^~ 2 as a function of —t for different values of y/s using the description according 
to Eq. (|15[) : blue solid curve: \/s = 3.686 GeV; green dashed curve: y/s = 10 GeV; the data correspond to the cross section 
measured at y/s = 3.686 GeV [la ]. 




FIG. 6. Kinematics of the pp — > nee process. 



Therefore, we choose the angles e + e - and $ e + e - as two independent kinematic variables describing the leptonic 
subprocess. The hadronic part of the amplitude depends on the cm. -energy y/s, the virtuality of the photon q 2 , and 
the Mandelstam variable t, which can be related to the pion-scattering angle 9 V . 
The differential cross section of the reaction is given by 

' W(0 e+e -, *»+.-), (22) 



dtdq 2 dn e+e - I6n 2 s(s- 4m 2 N ) (47r) 2 8g 2 3 
with the leptonic solid angle dfl e + e - . The e + e~ decay angular distribution W can be expressed as 

3 r 

W{0 e+e - , $ e+e - ) =— sin 2 6 e+e - poo + (1 + cos 2 6 e+e - )p u 

+ \/2 sin20 e + e - cos $ e + e -Re[pio] + sin 2 9 e + e - cos 2$ e+e - Re[pi_i] , (23) 
where the density matrix p is given by 



P\\> 



(M^,e;(X^))\ for A 7 , A; = 0,±1. (24) 



The expression presented in Eq. (|22p is model independent. However, in order to obtain numerical results we have 
to use a model to characterize the hadronic subprocess. We choose the Regge pole model which has been introduced 
m Section HU Since a virtual photon is produced, we have to specify the electromagnetic form factors parametrizing 
the 7* N N- vertices and 7* ./VA-vertices. For the electromagnetic form factors we consider a vector meson dominance 
(VMD) model, which assumes that the electromagnetic interaction is described through the exchange of the lowest 
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FIG. 7. Differential cross section da/dtdq 2 dQ e + e - of pp — > TT°e + e~ as a function of cos# e + e -. Red curves correspond to 
the TV-trajectory contribution; red solid curve: <& e + e - = 0; red dashed curve: $ e + e - = tt. Blue curves correspond to the 
(TV + A)-trajectory contribution including a reduction of the A pole residue (TF m 0.5) according to Eq. JT5J; blue solid curve: 
<J> e + e - = 0; blue dashed curve: $ e + e - = tt. 



lying vector mesons p, u, and <ft. In [2ll | a parametrization for both form factors F\ and F2 in the spacelike as well as in 
the timelike region is presented, which we will use for the purpose of the computation of the cross section. Eventually, 
the aim of the present work is to provide a further constraint on future extractions of timelike nucleon form factors. 

For the ma gne tic dipole form factor of the TV — > A transition, we use the results in the large N c limit, which can 
be written as j22j : 

g„tf) = iMM [F *tf)- F »tf) ]l (25) 

where F% (F£) is the Pauli form factors of the proton (neutron), for which we use the description of the VMD model, 
and Ky = K p — K„ = 3.70. 

The results of the differential cross section da/dtdq 2 dYL e + e - as a function of cos# e + e - are presented in Fig. [7]for 
several kinematical conditions. We display the iV-trajectory and (TV + A)-trajectory contributions as introduced in 
Eq. (fT5")) for the angles & e + e - = and $ e +e- = ^- The dependence on the angle $ e + e - appears as an asymmetric 
behavior of the cross section with respect to co§6 e + e -. For $ e + e - = it/2 the resulting cross section is symmetric, 
which can be derived from the general form of the decay angular distribution W, given by Eq. (|23[) . 

Using the Born model suggested in Q , one obtains a cross section, which is 1 to 4 orders of magnitudes larger than 
the results of the Regge pole model, depending on the variation of the kinematic parameters s, q 2 and 6^. 

The integrated cross section can be used to access the density matrix p. The cross section integrated over the 
azimuthal angle $ e + e - 

d<? f 2lT ,^ da 

» (26) 



dcos 6k dq 2 dcos e + e ~ Jo dtdq 2 dQ 
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FIG. 8. Cross sections denned by Eq. (|26[) (left) and Eq. (J27I) (right): Red (dashed) curve: iV-trajectory contribution; blue 
(solid) (JV + A)-trajectory contribution according to Eq. ()15p . 



is sensitive to poo and Pn and the cross section integrated over the polar angle 6 e + e - 

da f 1 da 

- — , — = / rfcosfl e+e - , (27) 

dcosbv dq z a3> e + e - V-i dtdq z dil e + e - 

can be analyzed in order to extract p\-\ from the $ dependence of the cross section. As selective predictions, we 
show the results of the cross sections (|26|) . and (|2~T)) in Fig. [5] for s = 5 GcV 2 , q 2 — 1 GeV 2 and 6^ = 7r/3 using both 
7V-tajectory and (N + A)-trajectory exchange, given by Eq. f|15|> . 



IV. CONCLUSIONS 



In this paper we studied the process pp — > 7r°e + e _ , giving a model independent expression of the cross section in 
terms of the lepton pair angular distribution and presenting results within a Reggc pole approach. We found that a 
model based on nucleon and A-Regge trajectory exchange provides a good description of the data for the process of 
real photoproduction pp — s- 7r° 7 in the energy range of s ~ 8.5 — 14GcV 2 . Applying this model to pp — ?> 7r°e + e _ allowed 
us to predict the angular dependence of the differential cross section, which can be used to extract the timclikc form 
factors in the unphysical region as well as their phases, in kinematics which will be accessible by the PANDA@FAIR 
experiment. 
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